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Abstract 

The effects of quantum lattice fluctuations on the Peierls transition and the opti- 
cal conductivity in the one-dimensional Holstein model of spinless fermions have been 
studied by developing an analytical approach, based on the unitary transformation 
method. We show that when the electron-phonon coupling constant decreases to a 
finite critical value the Peierls dimerization is destroyed by the quantum lattice fluctu- 
ations. The dimerization gap is much more reduced by the quantum lattice fluctuations 
than the phonon order parameter. The calculated optical conductivity does not have 
the inverse-square-root singularity but have a peak above the gap edge and there ex- 
ists a significant tail below the peak. The peak of optical-conductivity spectrum is not 
directly corresponding to the dimerized gap. Our results of the phase diagram and the 
spectral-weight function agree with those of the density matrix renormalization group 
and the exact diagonalization methods. 
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A great deal of quasi-one-dimensional materials, for example, the halogen-bridged mixed- 
valence transition-metal complexes, the conducting polymers, the organic and inorganic spin- 
Peierls systems, exhibit an instability against a periodic lattice distortion due to the Peierls 
dynamics. Among the models for one-dimensional systems the Holstein Hamiltonian|]l| is a 
typical electron-phonon coupling model studied by many previous authors. An interesting 
and still controversial problem is how the dimerized ground state is modified when quantum 
lattice fluctuations are taken into account. The quantum lattice fluctuations could have 
an important effect in most quasi-one- dimensional materials with a dimerized ground state 
because the lattice zero-point motion is often comparable to the amplitude of the Peierls 
distortion^]. The challenge of understanding the physics of quantum lattice fluctuations 
led to an intense study of the Holstein model. Generally speaking, the nonadiabatic effect 
suppresses the order parameters of the system 0]. As the optical-absorption concerned, the 
results of adiabatic approximations have inverse-square-root singularity at the gap edge. 
However, this approach is questionable and it has been shown that the quantum lattice 
fluctuations must be taken into account to satisfactorily describe some physical properties 
of quasi-one-dimensional systems By considering the nonadiabatic effect, the singularity 
may disappear |J. The influences of the phonon frequency on the optical-conductivity spec- 
trum in the range from cj q = to uj —>■ oo should be studied for understanding the physics 
of electron-phonon interactions in nonadiabatic case. 

When the quantum lattice fluctuations are taken into account the theoretical analysis 
becomes much more difficult. In the past several years, the Holstein Hamiltonian has been 
investigated by using various numerical approaches, such as Green's function Monte Carlo 
simulation H |7|], renormalization group analysis[§, [|, variational method of the squeezed- 
polaron wave-function [|lO|l , phenomenological approach^], exact diagonalizationfTT|, etc. Very 
recently, works of numerical approaches have been performed in relation to the Peierls transi- 
tion and the optical conductivity in the one-dimensional Holstein model of spinless fermions 
by using the density matrix renormalization group [I2| and the exact diagonalization[|T3[ 
methods. However, as was pointed out in Ref . [0 , because of the effects of limited system 
sizes in numerical approaches, the precise determination of the critical value in the small 
u regime is somewhat difficult and the precise extraction of the dimerized gap from op- 
tical conductivity data is prevented. An analytical study of the Holstein model will make 
it possible to have an insight into the intrinsical properties of the molecular crystal mate- 
rials. In a resent work, two of us[14| studied this model and investigated the dimerization 
order parameters and density of states in gapped phase and the velocity of charge excita- 
tions and Luttinger liquid stiffness constant in gapless phase. In this work we concentrate 
on the properties of the phase transition and the optical-responses of the system with the 
view of understanding the effects of quantum lattice fluctuations on the Peierls instability 
and the optical conductivity in the Holstein model. We will show that our results of the 
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phase diagram and the spectral-weight function agree surprisingly well with those of the 
density matrix renormalization group [12 and the exact diagonalizationfl~3|| methods, and in 



our theory the critical value can be determined precisely even for extremely small phonon 
frequency. The effects of quantum lattice fluctuations on the dimerized gap and on the or- 
der parameters are essentially different. The peak of optical-conductivity spectrum is not 
directly corresponding to the dimerization gap. 

The one-dimensional spinless Holstein model, in momentum space is 

H = Y, UJ V h \ h 1+Y, e k C W- -7^Y,9{ b q + b -q) C t+q C k, CO 



q k v N qik 

where e k = —2tcosk is the bare band structure, t the hopping integral, and N the total 
number of sites. c k and b q are the annihilation operators of electrons with momentum k and 



phonons with momentum q, respectively. The dispersionless phonon frequency u = J K/M 
and g is the electron-phonon coupling, K the elastic constant and M the mass of ions 
(throughout this paper we set % = ks = 1)- 

In order to take into account the electron-phonon correlation, the unitary transformation 



approach is used to treat ff||14j|, H = e s He s . After averaging the transformed Hamiltonian 



over the phonon vacuum state we get an effective Hamiltonian for the fermions 
#eff = \KNul + Y J Ev(k)c\c k - ]T A (k)(cl_ n c k + 4c fc _ ff ) 

Z k k>0 

-Jj E + ?. k )i 2 - W - <?> k')]cl +q c k cl_ q c k ,, (2) 

iV q,k,k' U 

where 

E {k) = e k - ±--£^6{k',k)6{k,k')(e k -e k ,), (3) 

A {k) = au [l-S{k-TT,k)}, (4) 

a = gy/2Muo, and S(k + q,k) = 1/(1 + \e k + q — ejfe|/k>o) is a function of the energies of the 
incoming and outgoing fermions in the electron-phonon scattering process. This effective 
Hamiltonian works well in the Uq = and ujq —>■ oo limits. 
u can be determined by the variational principle, 

U ° = ~Vn ^ 1 ~ 6 ^ k ^ 7T ' (f e | ( c I-7T C fc + cW-n) | fe) . (5) 

Here \fe) is the ground state of H e g. Thus, the total Hamiltonian can be written as H = 
H + Hi, where Hi includes the terms which are zero after being averaged over the phonon 
vacuum state, and 

#0 = + #eff. (6) 
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By means of the Green's function method to implement the perturbation treatment on the 
four-fermion term in Eq.([|), we get the renormalized band function and the gap function [jl4}| 

E k = E (k) — £{£(&', k)[2 — 8(k', k)] - S(k' - tt, k)[2 - 8(k' - ir, k)]} 

7V k'>0 



Afc = auo[c — d8(k — 7T,k)). (8) 



Where, the dimensionless coupling A = g 2 /2tu , W k = y E% + /S? k is the fermionic spectrum 
in the gapped state, and 

c = l + ^Y,m-7r,k)-V] , Mk) (9) 

N to auoJE^ + AKk)' 



d = i-^Et 1 -^-^^)] — / Ao(k) (io) 

N *>o aujE$(k) + Al(k) 



V = ]p E S(k,k + q)[2-S(k',k'-q)} (11) 

q,k,k' 

is the on-site interaction which should be subtracted because of the Pauli principle. 
The phonon-staggered ordering parameter is 

1 



m 



p N 



i 



a Eli- m 



These are basic equations in our theory. 

From Eq.(||), let Uo = 0, we get the self-consistent equation of phase transition points in 
the g 2 /u ~ u plane, 

1 4tA Y^h xfi , ^ c-dSjk - 7r,fc) 
N k>o \ h h\ 

If Uq — we have 5(k',k) = and c = 1, Eq.([5]) becomes the same as that in the 
adiabatic theory. In our theory 5(k — n, k) has a sharp peak at the Fermi point and, since 
1 — S(k — iT,k) = 4t| cos k\/(uo + 4t\ cos fc|), the logarithmic singularity in the integration of 
Eq.(^|) in the adiabatic case is removed as long as the ratio Wo/t is finite. Comparing Eq.flD 
with that in the adiabatic case, A = auo, we have the gap in the nonadiabatic case, 

A = A(7r/2) = auo[c — d\. (14) 

This is the true gap in the fermionic spectrum. 
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Fig.l shows the ground state phase diagram in the g/uj ~ t/uj plane. The solid line is 
our analytical result. For comparison, the results of previous authors are also shown. The line 
with circle and the line with square denote the results of the density matrix renormalization 
group (DMRG)[12[ and the two-cutoff renormalization group(TCRG) P3| , respectively. The 
line with triangle is the phase boundary of the variational Lanczos approach based on the 
inhomogeneous modified variational Lang-Firsov transformation (IMVLF)Jl3|. To check the 
consistency of the phase transition quantitatively with that of DMRG, some critical values 
g c are listed in Table |. One can see from both the figure and the table that our results agree 
surprisingly well with that of DMRG except for very large ujq (uj /t > 10). Furthermore, our 
theory can get the phase boundary, separating Luttinger liquid (LL) and insulation (CDW) 
phases, even in the very small ujq regime, which is theoretically and experimentally significant 
since, from the view point of experiment, for quite a lot of realistic cases the frequency of 
quantum phonon uq is small. It seems that in the DMRG and the finite-lattice Lanczos 
approach, because of the effects of limited system sizes, the precise determination of the 
critical value in the small ujo regime is somewhat difficult (JH| . The infinite system is never 
really gapless within the adiabatic approach, because the gap remains nonzero, although it 
becomes very small for weak electron-phonon coupling. On the contrary, in our theory, the 
logarithmic singularity J^ 2 dk/ cos A; in the integration of Eq.([5]) is removed by the factor 
1 — 5(k — 7r, k) and the critical value A c can be determined precisely even for extremely small 
phonon frequency. 

Inspired by the success of obtaining the phase diagram, we further investigate the dimer- 
ization gap and the optical responses. 

Fig. 2 shows the dimerization gap A/t = A(n/2)/t and the phonon order parameter 
arrip/t as functions of the phonon frequency uo/t for A = 0.81. It is the most notable that 
there is a discontinuous drop in the dimerization gap once the phonon frequency changes, no 
matter how small it is, from zero to finite, though at the adiabatic limit the dimerization gap 
A(n/2)/t = arrip/t. After the drop the dimerization gap and the phonon order parameter 
decrease as the phonon frequency increases. At the critical value u c , the dimerization gap and 
the phonon order parameter go to zero simultaneously and the system becomes gapless, which 
indicates that the quantum lattice fluctuations can destroy the dimerized Peierls state. The 
dimerization gap is much more reduced by the quantum lattice fluctuations than the phonon 
order parameter. The reason for the different behavior of the dimerization gap and the 
order parameter is that the former is the value of Eq. ffPf) at the Fermi point k = ir/2, where 
the quantum lattice fluctuations have the strongest effect, while the latter is the integral 
(see Eq.(JfJ)) over the all Brillouin zone and the effect of the quantum lattice fluctuations 
is gentled. The effects of quantum lattice fluctuations on the dimerization gap and on the 
phonon order parameter are essentially different, especially when uq is small. In the mean 
field (MF) approximation the Peierls distortion opens a gap 2Amf and Amf = ocm p . This 



5 



relation is sometimes assumed remains valid when quantum lattice fluctuations are taken 
into account. Our results indicate that this relation holds only in the adiabatic limit. 

The optical conductivity a {to) can be expressed by the retarded Green's function as 
follows: 



where K R is defined as 

K R (u) 



a{u) 



-/ / e 

-OQ 



2eonc 



ixuo 



lmK R (u), 



lUit 



dt{g\[j(0)j(t)-j(t)j(0)}\g). 



Here j is the current operator [16 



(15) 



(16) 



;i7) 



and j(t) = exp(iHt)j exp(— iHt) is the form of j in the Heisenberg representation. The 
unitary transformation of current operator is e s je~ s = j + [S,j] + ^[S, [S,j]} + 0(g 3 ). All 
terms of higher order than g 2 will be omitted in the following treatment. Because the 
averaging of Hi over the phonon vacuum state is zero, in the ground state at zero temperature 
Hi can be neglected. By using the approximately decoupling \g') « \g' Q ), the ground state 
of H , and H w H Q in the calculation 



{g\mm\9) 



(^|[e (S+fi) je-( s+ ^]e^[e( 5+ ^e-( s+ ^]e-^|^) 
(g' \[e s je- s ]e^Yje- S ]e-^ U), 



18) 



we can get 
K R (u) = 



+ 



fc>0 

J 2 



1 



J 2 

— y , 

N ~ \u - 2W k + i0~ 



1 



E 



N2 k>w>o u o V^-^o 



u + 2W k - i(F 
1 



sin k 



2^g 2 



W 2 



W k - Wy + iO" 
x[S 2 (k', k)(smk' - sin k) 2 (a k /3 k , + f3 k a kl ) 2 



to + u + W k + W k > - ^ 



+5 (k' - it, k) (sin k' + sin k) (a k a k > + /3 k f3 k 



(19) 



where a k 



1 + E k /W k )/2, and (3 k = J(l - E k /W k )/2. Thus, the optical conductivity 



crico 



2e ncJ 2 



+ 



- 2W k) sin 2 k 



k>0 

2er ) ncJ 2 



iV fc' u o 



5] ^5(cj - cu - Vr fc - WV)^', A;) (sin fc' - sin k) 2 (a k (3 k , + f3 k a k > 



uN 



k>0,k'>0 ^0 



+5 (k' — 7i, k) (sin k' + sin k) (a k a k / + (3 k f3 k 
and the u;— integrated spectral- weight function 



a{uj')duj' . 



(20) 



(21) 
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The optical conductivity for different phonon frequencies are shown in Fig. 3. The pa- 
rameter values used are: A = 1.0, with uo/t = 0.01,0.05, and 0.10. One can see that as 
u) increases, the optical-absorption spectrum broadens but the peak height decreases and 
moves to lower photon energy, and the spectral-weight increases as Uq increases. The inverse- 
square-root singularity at the gap edge in the adiabatic case|17]] disappears and there is a 
significant tail below the peak because of the nonadiabatic effect. We note that in our theory, 
in mathematical viewpoint, the difference between the ujq = and ujq > cases mainly comes 
from the functional form of the gap [see Eq.([8|)]. Comparing it with that in the adiabatic 
limit, one can see that the subgap states come from the quantum lattice fluctuations, i.e., 
the second term in the square bracket of Eq. (§) . 

The rescaled (^-integrated spectral- weight function S(u) / S m versus the photon energy uj/t 
relations of our result (solid line) and that of IMVLF (dashed line) are shown in Fig. 4, where 
5m = S{uj — > oo). The optical conductivity u(uj) (rescaled by S m ) of our result (solid line) is 
also shown. Because of the finite-size effects the optical conductivity of IMVLF is oscillatory 
and can not be compared with our result directly. We use the same parameter values as in 
Ref.0 (Fig.6(b)): g/u = 4.47, and u /t = 0.1. The true gap 2A(vr/2)/t = 1.04 obtained 
from Eq . (|T4]) is marked by the arrow. One can see that the peak of optical-conductivity 
spectrum is not directly corresponding to the dimerized gap. The energy gap is small than 
the activation energy of the optical-conductivity. Our result shows clearly both the position 
and the peak of optical conductivity and the spectral-weight agrees with that of IMVLF. 
While in the exact diagonalization method the finite-size effects prevent a precise extraction 
of the CDW gap from the optical exact diagonalization data[[L3|. 

In conclusion, the effects of quantum lattice fluctuations on the optical-conductivity 
spectrum and the ground state phase diagram of the one-dimensional Holstein model of 
spinless fermions have been studied by developing an analytical approach. We show that 
when the electron-phonon coupling constant decreases the dimerization gap decreases, and at 
a finite critical value the Peierls dimerization is destroyed by the quantum lattice fluctuations. 
The critical value of electron-phonon coupling can be determined precisely even for very 
small phonon frequency. The dimerization gap is much more reduced by the quantum 
lattice fluctuations than the order parameter. The calculated optical conductivity does 
not have the inverse-square-root singularity but have a peak above the gap edge and there 
exists a significant tail below the peak. In nonadiabatic case the dimerization gap is small 
than what the peak position of the optical conductivity is corresponding to. Our results of 
the phase diagram and the spectral-weight function agree with those of the density matrix 
renormalization group and the exact diagonalization methods. 
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Figure Caption 



Fig.l The ground state phase diagram in the g/w ~ t/u plane. The solid line is our 
analytical result. The line with open circle and the line with open square denote the results 
of the density matrix renormalization group (DMRG) and the two-cutoff renormalization 
group(TCRG) methods, respectively. The line with open triangle is the phase boundary of 
the variational Lanczos approach (IMVLF). 

Fig. 2 The dimerization gap A/t = A(n/2)/t and the order parameter am p /t as func- 
tions of the phonon frequency Do/t in the case of A = 0.81. 

Fig. 3 The optical conductivity in the case of A = 1.0 for different phonon frequencies 
wo/t = 0.01,0.05, and 0.10. 

Fig. 4 The rescaled optical conductivity (t(u) and tu-integrated spectral- weight function 
S(u)/S m versus the photon energy ui/t relations of our results and that of IMVLF, where 

= S{uj — > oo). The parameter values used are same as in Ref.fl4|1 (Fig.6(b)): g/uio — 4.47, 
and uio/t = 0.1. The dashed line is the result of IMVLF and the solid lines are our results. 
The arrow marks the true gap 2A(ir/2)/t = 1.04. 
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Table Caption 



Table I: Critical point g c . g* is the value of g determined by letting the stiffness constant 
K p = \. 
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